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ABSTRACT 

The X-radiation of large-scale extragalactic jets poses serious challenge for con- 
ventional electron-synchrotron or inverse Compton models suggested to explain the 
overall nonthermal emission of the resolved knots and hot spots. In this paper I pro- 
pose an alternative mechanism for X-ray emission - synchrotron radiation by extremely 
high energy protons - and discuss implications of this model for the extended jet fea- 
tures resolved by Chandra in several prominent radiogalaxies and AGN like Pictor A, 
3C 120, PKS 0637-752, and 3C 273. I show that if protons are indeed accelerated to 
energies E-p > 10^* eV, it is possible to construct a realistic model which allows an ef- 
fective cooling of protons via synchrotron radiation on quite "comfortable" timescales 
of about 10"^ - 10^ yr, i.e. on timescales which provide effective propagation of pro- 
tons over the jet structures on kpc scales. This explains quite naturally the diffuse 
character of the observed X-ray emission, as well as the broad range of spectral X-ray 
indices observed from different objects. Yet, as long as the proton synchrotron cooling 
time is comparable with both the particle escape time and the age of the jet, the 
proton-synchrotron model offers an adequate radiation efficiency. The model requires 
relatively large magnetic field of about 1 mG, and proton acceleration rates ranging 
from Lp ~ 10'*'^ to 10''^^ erg/s. These numbers could be reduced significantly if the 
jet structures are moving relativistically towards the observer. I discuss also possible 
contributions of synchrotron radiation by secondary electrons produced at interactions 
of relatively low energy [Ep < 10^^ eV) protons with the compressed gas in the jet 
structures. This is an interesting possibility which however requires a very large prod- 
uct of the ambient gas density and total amount of accelerated protons. Therefore it 
could be treated as a viable working hypothesis only if one can reduce the intrinsic X- 
ray luminosities assuming that the regions of nonthermal radiation are relativistically 
moving condensations with Doppler factors 5j ^ 1. The kpc scales of knots and hot 
spots are not sufficient for effective confinement of > 10^^ eV protons. This suppresses 
the synchrotron radiation by secondary electrons produced at p7 interactions. At the 
same time the runaway protons interacting with 2.7 K cosmic microwave background 
radiation, initiate non-negligible diffuse X- and 7-ray emission in the surrounding clus- 
ter environments. I discuss the spectral and angular characteristics of this radiation 
which essentially depend on the strength of the ambient magnetic field. 

Key words: radiation mechanisms: non-thermal - galaxies: jets - galax- 
ies:active:individual: Pictor A, 3C 120, 3C 273, PKS 0637-752 



1 INTRODUCTION 

The recent Chandra observations revealed various bright X- 
ray features like knots and hot spots in large-scale extra- 
galactic jets (for a review see e.g. Harris 2001). The undis- 
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puted synchrotron origin of the radio emission of the re- 
solved jet structures in several powerful AGN and radio- 
galaxies implies presence of relativistic electrons with a typ- 

o 1/2 1/2 

ical Lorentz factor 7e — 10 ^^ghz' ^h^re SmG = 

-B/10~"^G is the magnetic field and vqhz = i^/1 GHz is 
the frequency of synchrotron radiation. The inverse Comp- 
ton scattering of same electrons leads to the second compo- 
nent of nonthermal radiation. The 2.7 K cosmic microwave 
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background radiation (CMBR) provides a universal photon 
target field for the inverse Compton emission. The electrons 
boost the energy of the seed photons to _E ~ AhTj^; this 
gives a simple relation between the energies of inverse Comp- 
ton and synchrotron photons produced by the same electron: 
iJ ~ 1 i^GHzS~Q keV. Thus, for the characteristic magnetic 
field, which in the knots and hot spots is believed to be 
between 0.1 mG and 1 mG (e.g. Harris 2001), the inverse 
Compton radiation of radioelectrons appears in the Chandra 
energy domain. Yet, the energy density of CMBR is too low 
to provide significant X-radiation. Indeed, the inverse Comp- 
ton luminosity constitutes only a small part of the emission 
of radioelectrons, Lx ~ 10~^-B~qLh,, unless we assume that 
the radio emitting regions are relativistically moving con- 
densations with the Lorentz factor of bulk motion F > 10, 
thus the energy density of CMBR "seen" by electrons in the 
jet's frame is enhanced by a factor of (Tavecchio et al. 
2000; Celotti et al. 2001). This hypothesis automatically im- 
plies that the jet features are moving towards the observer, 
otherwise we would face unacceptably large requirements to 
the energy budget of the central source. 

The emissivity of inverse Compton radiation could be 
significantly enhanced also due to additional seed photon 
sources, in particular due to the jet's own synchrotron 
radiation. The density of this radiation in a knot at 
a distance d from the observer is estimated as uir — 
{d/Rfc-^fR-o =i 0.9 f-i2e-^ eV/cm^ where /_i2 = 
/r-o/10~^^ erg/cm'^s is the observed radio-to-optical flux 
normalized to lO"^'^ erg/cm^s, R is the source radius and 
6 — R/d is the source angular size in arcseconds. For ex- 
ample, the fluxes of 1 to 1000 GHz radio synchrotron emis- 
sion of two distinct hot spots ( "A" and "D" ) of the radio- 
galaxy Cygnus A, both of ~ 0."5 angular size, are very high, 
/ ~ 10"" erg/cm^s (Wilson et al. 2000), thus the density of 
synchrotron radiation exceeds the density of CMBR by two 
orders of magnitude. Therefore, the X-ray emission of these 
hot spots can be readily explained by the synchrotron-self- 
Compton (SSC) model, as it was predicted by Harris et al. 
(1994), assuming B ~ 0.2 mG, and taking into account that 
the energy flux at radio frequencies exceeds the X-ray flux 
by 1.5 orders of magnitude. 

This is not the case, however, for the majority of jet 
features resolved by Chandra for which we observe, in fact, 
just an opposite picture. The X-ray fluxes, for example, from 
the jets in Pictor A (Wilson et al. 2001), 3C 120 (Harris et 
al. 1999), 3C 273 (Roser et al. 2000; Sambruna et al. 2001; 
MarshaU et al. 2001), and PKS 0637-752 (Schwartz et al. 
2000; Chartas et al. 2000) are at least factor of 10 larger 
than the radio and optical fluxes. 

For X-ray emission of many resolved extragalactic knots 
and hot spots the electron synchrotron radiation is consid- 
ered as 'the process of choice' (Harris 2001). The synchrotron 
radiation is indeed an extremely effective mechanism con- 
verting the kinetic energy of relativistic electrons into X-rays 
with an almost 100 per cent efficiency, provided, of course, 
that the electrons are accelerated to multi-TeV energies. The 
single power-law spectrum smoothly connecting the radio, 
optical and X-ray fluxes with spectral index a ~ 0.76 ± 0.02 
observed from the knot "Al" in the jet of 3C 273 (Roser 
et al. 2000; MarshaU et al. 2000) formally could be brought 
as an additional argument for the synchrotron origin of X- 
rays. Note, however, that the observed spectral steepening 



(or cutoff) at optical frequencies (but with yet hard X-ray 
spectra) in the so-called "problem sources" (Harris 2001) 
like the knot D/H3 in the same 3C 273 (Roser et al. 2000, 
Marshall et al. 2000) , as well as the knots in 3C 120 (Harris 
et al. 1999), PKS 0637-752 (Schwartz et al 2000) and in the 
hot spot of Pictor A (Wilson et al. 2000), tell us that we 
deal with a rather complex picture compared, in particular, 
with the simple single-power-law synchrotron model. 

Generally, the high radiation efficiency is treated as a 
key component of any successful model/mechanism of lumi- 
nous nonthermal sources. However, for the large-scale ex- 
tragalactic jets this highly desired feature ironically leads to 
certain difficulties, as long as it concerns the synchrotron 
X-ray emission of electrons. The latter, in fact, seems to 
be "over-efficient". This, at first glance paradoxical state- 
ment, has a simple explanation. The cooling time of elec- 
trons responsible for radiation of an X-ray photon of energy 

yr is small (for any rea- 
sonable magnetic field B > 0.01 mG), even compared with 
the minimum available time - the light travel time across 
the source, R/c ~ 3 x 10^ yr. This simple estimate has the 
following interesting implications. 

(a) The acceleration of electrons should take place 
throughout entire volume of the hot spot, otherwise the 
short propagation lengths of electrons would not allow for- 
mation of diffuse X-ray emitting regions on kpc scales as it 
is observed by Chandra. However, operation of huge quasi 
homogeneous 3-dimensional accelerators of electrons with 
linear size of 1 kpc or so, seems to be a serious theoretical 
challenge. Possible alternatives could be (1) a hypothesis of 
superposition of many but compact accelerators within the 
observed knots, or (2) a scenario when the cloud of TeV 
electrons, accelerated in a relatively compact region, expand 
with the speed of light, e.g. in the form of a (quasi) sym- 
metric relativistic "hot" windQ 

(b) We should expect quick formation of the spectrum of 
cooled electrons with power-law index Fc = Fo 4- 1, where 
Fo is the acceleration spectrum index. The power-law index 
of radiating electrons Fc is found from the spectral index 
of synchrotron radiation a {Su oc as Fo = 2q -f 1, 

therefore Fo = 2a. For some extragalactic jets, in particular 
for the knot "Al" of 3C 273 with a power-law spectral in- 
dex in the X-ray band Qx ~ 0.6 (Marshall et al. 2001), this 
would imply an extremely hard acceleration spectrum with 
Fo = 1.2 ± 0.1 (!) which would challenge any (in particular, 
shock) acceleration scenario. Moreover, the reported flux up- 
per limit at the optical band (< 0.18/iJy at 6.8 x 10^** Hz) 
compared with the X-ray flux (O.OlS^Jy at 2 keV) from the 
so-called 25" knot of 3C 120 (Harris et al. 1999) implies that 
the X-ray to optical spectrum is flatter than Si, oc u'^'^^, 
and therefore the spectrum of radiating electrons is flatter 
than N{Ec) oc E~^''' . Harris (2001) has interpreted this re- 
sult as a deviation from the canonical shock acceleration 
scenario but still in accordance with oblique shock acceler- 
ation models (e.g. Kirk 1997). However, the E^^ type dif- 
ferential electron spectrum is, in fact, the flattest possible 



t The "hot wind" means that the electrons in the frame of 
the wind have a broad relativistic momentum distribution which 
would enable their synchrotron radiation. 
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(independent of the acceleration model) cooled steady-state 
spectrum of electrons, as long as the synchrotron (and/or 
inverse Compton) energy losses play dominant role in its 
formation. Therefore, in order to explain the observed spec- 
trum of synchrotron radiation with Qo-x < 0.5 wc must 
keep the acceleration spectrum Q{Ee) unchanged, i.e. avoid 
radiative losses. Formally this could be possible within the 
above mentioned hypothesis of fast, with speed- of -light es- 
cape of electrons in a form of above mentioned "hot" rela- 
tivistic wind. This is, however, an ad hoc assumption which 
requires thorough theoretical and observational inspections. 

In this paper we suggest a different approach for in- 
terpretation of the diffxisc nonthermal emission from largo 
scale extragalactic jets at X-rays and, perhaps, also at lower 
frequencies. Namely, we assume that this radiation has a 
synchrotron origin, but is produced by very high energy 
protons. We will show that adopting magnetic fields in jet 
structures somewhat larger (by a factor of 10 or so) than 
the field assumed in the standard electron synchrotron or 
inverse Compton models, and speculating that the protons 
are accelerated in the jet structures to energies Ep ~ 10^* eV 
or more, we can construct a model which allows an efficient 
cooling of protons via synchrotron radiation on very "com- 
fortable" timescales, namely on timescalos comparable with 
the (diffusive) escape time of protons and/or the age of the 
jet. This provides effective propagation of protons over the 
entire jet structures on kpc scales, and thus can naturally 
explain the diffuse character of X-ray emission, as well as 
the broad range of spectral indices of X-rays observed from 
different objects. Yet, as long as the proton synchrotron cool- 
ing time, the escape time of particles and the age of the jet 
are of the same order of magnitude, the proton synchrotron 
model offers a quite high (from 10 to almost 100 per cent) ef- 
ficiency of transformation of the kinetic energy of protons to 
hard synchrotron X-ray emission. This makes, in our view, 
the proton synchrotron radiation an attractive and viable 
mechanism for interpretation of nonthermal X-ray emission 
from large-scale extragalactic jets. 



2 EXTREMELY HIGH ENERGY PROTONS IN 
THE AGN JETS 

The jets of powerful radiogalaxies and AGN are one of a few 
potential sites in the Universe where protons could be accel- 
erated to highest observed energies of about 10^" oV (see e.g. 
Hillas 1984; Cesarsky 1992; Rachen & Biermann 1993; Henri 
et al. 1999). Dissipation of bulk kinetic energy and/or the 
Poynting flux of jets results in, most likely through strong 
terminal shocks, generation/amplification of magnetic fields, 
heating of the ambient plasma, and acceleration of particles 
in the knots and hot spots. More specifically, Rachen and 
Biermann (1993) have demonstrated that mildly- relativistic 
jot terminal shocks in hot spots of FRII radiogalaxies with 
typical magnetic field of about 0.5 mG and size ~ 1 kpc 
are indeed able to accelerate protons up to lO'^" cV. Os- 
trowski (1998) has shown that very effective acceleration of 
extremely high energy protons (with a rate of ~ 10 per cent 
of the maximum possible acceleration rate, t^^c ~ c/rg) can 
take place also at the jet shear boundary layers. 



2.1 Secondary electrons 

The undisputed synchrotron origin of nonthermal radiation 
of many knots and hot spots extending to (at least) optical 
wavelengths is treated as a clear evidence for acceleration of 
electrons to TeV energies. Besides, relativistic electrons are 
produced also in interactions of high energy protons with 
ambient gas and photon fields. Biermann and Strittmatter 
(1987) has proposed that the secondary component of elec- 
trons produced by extremely high energy protons interacting 
with local photon fields, may play an important role in for- 
mation of the nonthermal radiation of radio jots. This idea 
has been further developed within the PIC (Proton Induced 
Cascade) model by Mannheim et al. (1991) for interpreta- 
tion of radiation of some prominent radio jets features like 
the west hot spot in Pictor A and the knot Al in 3C 273. 
The intrinsic feature of this model is that it produces a dis- 
tinct maximum in the spectral energy distribution (SED) 
vSi, at MeV energies, and a standard hard X-ray spectrum 
with a spectral index ~ 0.5 at keV energies. The recent 
Chandra observations of Pictor A (Wilson et al. 2001) give 
a flux which at 1 keV appears almost 3 orders of magni- 
tude below the PIC-model predictions. The predictions of 
this model for the knot Al of 3C 273 also fall well below the 
reported Chandra flux at 1 keV (Sambruna et al. 2001; Mar- 
shall et al. 2001). Although it is formally possible to increase 
the X-ray fluxes, assuming significantly larger power in ac- 
celerated protons or speculating with denser target photon 
fields, this would however face serious problems with the re- 
quired nonthermal energy budget. Moreover, the flat X-ray 
spectrum of the hot spot of Pictor A with a spectral in- 
dex a = 1.07 ± 0.11 (Wilson et al. 2001) excludes the PIC 
model which predicts much harder X-ray spectrum. The X- 
ray spectral index of the knot Al in 3C 273 Ox = 0.60 ±0.05 
is closer to the PIC model predictions, but in this source the 
efficiency of the process is too low to explain the detected 
absolute X-ray flux. 

For a broad, e.g. power-law spectrum of target photons, 
the 'photo-meson' cooling time of protons can be estimated 
as tpj = [c < af > n(f*)c*]"-', where < af >~ 10"^* cm^ 
is the photo-meson production cross-section weighted by in- 
elasticity at the photon energy ~ 300 GeV (see e.g. Miicke 
et al. 1999) in the proton rest frame and e* — O.OSE^q eV; 
Eig = Ep/W^^ eV is the energy of the proton in units 
of lO"'^ eV. The fluxes of jet features at most relevant for 
photo-meson production frequency band 10^^- lO^'' Hz typi- 
cally are poorly known. Therefore it is convenient to present 
the low-frequency photon flux density in the following form 
5*1^ = SqVq^^, where is the flux density at 1 GHz in units 
of Jy. This gives 

tp^ ~ lOVaSo-'SiT^^ yr, (1) 

where Cc « 0.42; 3.9 and 36 for a = 0.5, 0.75 and 1.0, 

respectively; 8 = R/d is the angular radius of the source 
in arcseconds. For example, in the Al knot of 3C 273 with 
6 ~ 0.5 — 1 arcsecond, So — 0.05Jy and a ~ 0.75 (Marshall 
et al. 2000) , the characteristic photo- meson production time 
appears more than 2 x 10® yr even for 10^° eV protons. This 
is too large compared with the escape time of highest en- 
ergy protons from the X-ray production region. Indeed, the 
particle escape cannot be longer than the time determined 
by diffusion in the Bohm limit, iesc /2D with diffusion 
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coefficient D(E) = rjr^c/S, where rg = Ep/eB is the the 
gyroradius, and 77 > 1 is the so-called gyro-factor (in the 
Bohm limit rj = 1). Thus 

te.c ^ 4.2 X 10'7)-'B„,Gi?kpc-BiV yr- (2) 
Thus, representing the source radius in the form R = 9 ■ d '2:^ 
AM{d/l Gpc) kpc, for 3C 273 {z = 0.158) we findQ 

tescAp-, =i 8 X 1Q-\B^gE-°-^\ (3) 

This ratio implies a very limited efficiency of transforma- 
tion of the kinetic energy of protons to nonthermal radia- 
tion. Note that both characteristic times, tosc and tp-^ are 
proportional to R , consequently the ratio tosc/tp7 does not 
depend on 6, thus we cannot increase the photo-meson pro- 
duction efficiency assuming that the acceleration and radia- 
tion take place in smaller (but numerous) regions inside the 
knot. The ratio slightly depends also on the energy of pro- 
tons. Although the photo-meson production time decreases 
with energy, the escape time decreases even faster, thus the 
efficiency of the process (tasc/tp-y) becomes, in fact, less at 
higher energies. The direct (Bethe-Heitler) production of 
{e'^,e~) pairs cannot significantly enhance the efficiency of 
P7 interactions; although this process requires lower energy 
protons which are confined more effectively, for a broad- 
band target radiation this process typically is slower than 
the photo-meson production. The fluxes of mm radiation in 
some hot spots, e.g. in the radiogalaxies Pictor A (Wilson et 
al. 2001) and Cygnus A (Wilson et al. 2001) are 2 or 3 orders 
of magnitude higher than in the knot Al of 3C 273, there- 
fore the efficiency of p7 interactions in these objects could 
approach to a few per cent, provided that particle escape 
takes place in the Bohm regime. Even so, the PIC remains a 
rather ineffective scenario for production of X-rays, because 
it allows less than 1 per cent of the overall luminosity of 
nonthermal radiation to be released at keV energies. 

The low-frequency radiation converts, through photon- 
photon interactions, a fraction of high energy 7-rays, before 
they escape the knots, into {e'^,e~) pairs. Thus, these in- 
teractions could initiate, in principle, electromagnetic cas- 
cades in the jet. However, it is easy to show that for most of 
the X-ray features detected by Chandra, this process could 
be neglected. Indeed, for the power-law spectrum of low- 
frequency radiation with above defined normalization, the 
optical depth which characterizes the efficiency of the cas- 
cade development at 7-ray energy E is estimated as 

T-y-, ^ T^Soe-^R^p.{E/l TeV)" , (4) 

with Tc, ~ 4 X 10"^; 1.8 x 10"'' and 8.2 x 10"'' for a =0.5, 0.75 
and 1.0, respectively. For example, in the knot A of 3C 273 
with ~ 0.5 arcsec, a ~ 0.75 and So — 0.05 Jy, the optical 
depth exceeds 1 only aX E > 10^^ eV. Thus, the cascade 
could be initiated only at the presence of very high energy 
7-rays. Even so, for the typical magnetic fields B > 0.1 
mC in the knots and hot spots, the Compton cooling of the 
secondary electrons is a much slower process compared with 
the synchrotron losses (wb 3> Wr), thus in these objects the 
cascade development is strongly suppressed. 

Secondary electrons can be produced also at interac- 
tions of relativistic protons with the ambient gas. The char- 

■t Hereafter for the Hubble constant we adopt Ho = 60 km/s Mpc. 



acteristic cooling time of this process in the medium with 
gas number density no is almost independent of the proton 
energy 

tpp = (nocTpp/c)"^ ~ 1.7 X 10*n"^ yr . (5) 

Because the X-ray synchrotron radiation in the jet could be 
produced by < lOTeV electrons, the pp interactions require 
relatively low energy protons, Ep < 100 TeV the escape 
time of which may exceed the jet age. On the other hand 
the pp cooling time of protons is quite large, therefore this 
attractive mechanism, which produces relativistic electrons 
throughout the entire jet structure, may have an impact on 
the overall nonthermal radiation, provided that the ambient 
gas density in the knots and hot spots exceeds 0.1 — 1 cm"''. 

2.2 Synchrotron radiation of protons 

During the jet lifetime most of the energetic particles es- 
cape the knots. In particular, in the Bohm regime, when 
the particles tend to drift away most slowly {rj = 1), the 
break in the initial proton spectrum takes place at Ep ~ 
4.2 X 10"BmGi?kpc(At/10^ yr)"^ eV. Since the observations 
limit the size of a typical knot or a hot spots by i? < few kpc, 
the only possibility to move the break point to higher ener- 
gies is the increase of the magnetic field, B > 1 mC. How- 
ever, at such large magnetic fields the proton synchrotron 
radiation (see e.g. Aharonian 2000) becomes an additional, 
or even more important limiting factor with characteristic 
cooling time 

feynhc ^ 1.4 X lO' B-^E^g' yT , (6) 

Remarkably, the radiation by highest energy protons is re- 
leased in the X-ray domain: 

hryn, = 0.29/i!/c ^ 2.5 x 10^BmG-B?9 keV. (7) 

Correspondingly, the characteristic time of radiation of a 
proton-synchrotron photon of energy E is 

t{E) ~ 2.2 X 10'B-^\E/1 keV)"^/2 yr . (8) 

Taking into account that 1 keV photons are produced by 
protons with energy Eig ~ 0.063 B^^ , the time of radi- 
ation of synchrotron X-ray photons becomes less than the 
time of particle escape in the Bohm regime if Bf^^QR^p^ > 30. 
For a typical size of jet knots of about 1 kpc, this would 
require B > 3 x 10^'^ G and correspondingly the total en- 
ergy contained in the the magnetic field Wb = B^R^/Q > 
4 X 10^* erg. Interestingly, in such an extreme regime, when 
almost the whole kinetic energy of accelerated protons is 
released in synchrotron X-rays for a typical time of about 
< 4 X 10^ yr, the luminosity of synchrotron radiation at 1 
keV is expected at the level Lx — 3 x 10*^k erg/s, assum- 
ing equipartition between the magnetic field and accelerated 
protons, Wp — Wb - The parameter k is the fraction of the ki- 
netic energy in protons responsible for synchrotron X-rays, 
-Ep > 10^^ eV, to the total energy. For a Ep'^ type pro- 
ton spectrum extending beyond lO'^* eV, k ~ 0.1. Note that 
in the regime dominated by proton synchrotron losses, the 
total luminosity of the proton-synchrotron radiation does 
not depend on the magnetic field, and is determined sim- 
ply by the acceleration power of highest energy protons, 
Lx — Lp{> Ep). 



© 0000 RAS, MNRAS 000, 000-000 



Proton synchrotron radiation 5 



Below we present detailed numerical calculations for 
several prominent extragalactic jet features. We assume that 
during the lifetime of the jet At, protons are injected into 
a homogeneous and spherically symmetric region, with ac- 
celeration spectrum Qp{Ep) = Qo-Ep ^ exp {—Ep/Eo) and 
(quasi) continuous constant rate Lp = J Qp{Ep)EpdEp. The 
kinetic equation for energy distributions of relativistic parti- 
cles which takes into account the radiative and escape losses, 
and the time- dependent solutions to this equation, are de- 
scribed in the papers by Atoyan & Aharonian (1999) and 
Aharonian & Atoyan (1999). In this study the interactions 
of protons with the magnetic field {synchrotron radiation), 
photon fields {photo-meson- and Bethe-Heitler pair produc- 
tion) and the ambient gas {pp n —> e^e~7), as well as the 
synchrotron radiation and inverse Compton scattering of sec- 
ondary electrons are included in calculations. The fast syn- 
chrotron cooling of secondary electrons in relatively strong 
magnetic fields allow us to ignore the cascade processes (ini- 
tiated by 7-7 interactions) without any significant impact 
on accuracy of calculations. 



3 JET FEATURES IN PICTOR A, PKS 
0637-752, AND 3C 120 

The luminosities of proton-synchrotron radiation calculated 
for 3 prominent jet features - the western hot spot in Pictor 
A, the radio knot at 25" distance from the nucleus of the 
Seyfert galaxy 3C 120, and the bright X-ray region in the 
inner western jet of PKS 0637-752 (at ~ 8" from the core) 
- are shown in Fig. 1. The shape of the resulting spectrum 
of protons primarily depends on the index of acceleration 
spectrum Fo, the magnetic field B and the size R. The total 
power in the accelerated protons Lp determines the absolute 
X-ray flux. For clearness, in Fig. 1 we present the results of 
calculations which take into account interactions of protons 
with universal 2.7 K CMBR, but ignore the interactions with 
the ambient gas and local photon fields. 

3.1 Pictor A 

The western X-ray hot spot at 4'. 2 from the nucleus of 
this nearby (2 — 0.035) powerful radiogalaxy coincides with 
the radio jet and is laterally extended by ~ 2 kpc (Wil- 
son et al. 2001). The X-ray emission of the hot spot is 
well described by a flat power-law spectrum with a pho- 
ton index a + 1 — 2.07 ± 0.11. This implies that within 
the proton-synchrotron model the spectral index of radiat- 
ing protons should be close to 3. Thus, for the shock ac- 
celeration spectrum with canonical power-law index Fq ~ 2 
we must allow an increase of the spectral slope by 1. Such 
a steepening can be naturally provided by an intensive 
synchrotron cooling and/or by fast energy-dependent es- 
cape of protons from the source with Tcsc{E) oc 1/E. In 
the synchrotron-loss dominated regime the corresponding 
break in the proton-synchrotron spectrum (the point where 
Si, oc is changed to oc v^^) takes place at 

{hu)t ~ 5 B-^At^' keV , (9) 

where Atg = At/10* yr is the duration of the particle ac- 
celeration in units of 10* yr. For a typical lifetime of the jet 
(the period of activity of the central engine) At < 10* yr, 
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Figure 1. Proton synchrotron luminosities of jet features in Pic- 
tor A (western radio hot spot), 3C 120 (25" radio knot), and 
PKS 0637-752 (the outer jet component). The model parameters 
are described in the text. The optical and X-ray data are from 
Wilson et al. (2001) for Pictor A; Harris et al. (1999) for 3C 120; 
Schwartz et al. (2000) for PKS 0637-752. 

the flat {vSv = const) spectrum of the hot spot above 1 
keV would require magnetic fleld larger than 1 mG. In the 
regime dominated by particle escape, the break in the spec- 
trum depends on the propagation character of protons. In 
particular in the regime of Bohm diffusion the break in the 
proton-synchrotron spectrum appears at 

{hu)b ~ 0.004 Bl,GAt^^Rtp^ri~^ keV . (10) 

Note that the characteristic times of synchrotron cool- 
ing and the particle escape in the Bohm regime have the 
same oc 1/E dependence, therefore the dominance of syn- 
chrotron losses does not depend on the speciflc energy inter- 
val, and is determined by the following simple condition: 

bIbRIpcV'' > 30 . (11) 

In Fig. 1 we show two "proton-synchrotron" fits to the 
spectrum of the hot spot of Pictor A calculated (a) for a 
regime dominated by synchrotron losses with two key model 
parameters B =3 mG and R — 2 kpc (curve 1), and (b) for 
a regime dominated by particle escape with B = 1 mG and 
R = 1 kpc (curve 2). In both cases the following additional 
model parameters were assumed: gyro-factor r) = 1 (Bohm 
diffusion); power-law index Fo = 2 and maximum energy of 
accelerated protons Eo = 10^*^ eV. And finally, to match the 
absolute X-ray fluxes we assume that the proton acceleration 
takes place in a (quasi) continuous regime during last At = 
10* yrs with rates Lp = 4.5 x 10''^ and 1.35 x lO^'^erg/s for 
the cases (a) and (6), respectively. 

Although both spectra satisfactorily flt the reported X- 
ray data above 0.5 keV|^ the case (a) is more attractive 

§ Below 0.5 keV the spectrum corresponding to the case (a) 
passes below the range of extrapolated fluxes given by Wilson 
et al. (2001). However better agreement at these energies could 
be easily achieved assuming a slightly larger magnetic field. 
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Figure 2. Total energy of relativistic protons confined in X-ray 

emitting regions of radio jets: the western iiot spot in Pictor A 
(eurves 1 and 2), 25" knot of 3C 120 (curves 3 and 4), and tiie 
outer jet component in PKS 0637-752 (curves 5 and 6) . 

on energetic grounds - it requires ~ 30 times less power in 
accelerated protons. Note that in the regime dominated by 
synchrotron losses, the efficiency of X-ray production is close 
to 100 per cent, the luminosity being almost independent 
of the magnetic field, Lx — i ~ Lp. In the particle-escape 
dominated regime, -Lx-7 oc B^. 

The spectra of relativistic protons Wp{Ep) = EpN{Ep) 
trapped in the hot spot arc shown in Fig. 2. Like Fig. 1, 
the curves 1 and 2 correspond to (a) and (b) cases, respec- 
tively. The spectra are similar with an original (accelera- 
tion) shape oc Ep^ below the break energy at which the 
escape and synchrotron cooling times exceed the age of the 
source, and oc Ep^ at higher energies. At the same time, 
while in the case (a) the total energy in protons is rather 
modest, wi*°*^ = / W{Ep)AEp ~ 1.1 x 10^^ erg, in the 
case (6) it is significantly higher, Wj$'°'^ ~ 2.9 x 10''" erg. 
Thus, while in the case (o) the conditions are close to 
the equipartition between the protons and magnetic filed 
(wp ~ 1.3 X 10"'^ erg/cm^, wb ^ 3.6 x 10"'' erg/cm^s), 
in the case (6) the pressure of protons exceeds the B-field 
pressure by 3 orders of magnitude. 

Finally we note that the proton synchrotron radiation 
cannot be responsible for the observed radio and optical 
fluxes. The two optical/UV points (triangles) in Fig. 1 ob- 
tained from the HST observations of the hot spot (Wilson et 
al. 2001) confirm the break in the spectrum at ~ 10^* Hz 
reported earlier by Meisenheimer et al. (1997). Both the 
slope of the spectrum connecting the radio and optical 
fluxes, as well as the spectral break above 10^* Hz can be 
naturally explained by synchrotron radiation of electrons 
(Meisenheimer et al. 1997). We note, however, that this ra- 
diation component is produced, most probably, in region(s) 
separated from the X-ray production region, both in space 
and time. Indeed, the break in the spectrum of electron syn- 
chrotron radiation in the regime when the latter dominates 
in the formation of the spectrum of parent electrons, ap- 
pears at Vf, ~ 60 B~QA/;g ^ Hz. Thus, the synchrotron-loss 
origin of the break at 10^* Hz would require, for any rea- 
sonable magnetic field B > 0.01 niG, that the acceleration 
of this electron population took place recently, namely since 
last 10® years or less. Such a conclusion could be avoided if 
we assume that the electrons escape the source so fast, that 



the synchrotron losses become negligible. The comparison of 
the particle escape time 

tesc ^ 3.2 X lO^iikpc {c/v) yr (12) 

with the radiation time of a synchrotron photon of frequency 

p. 

til.) ^ 75B;Jf\u/lQ'^ Hz)-^/Vr (13) 

lead to the conclusion that even for the escape of electrons 
with the speed of light, the magnetic field should be less 
than 0.1 mG. 

The calculations presented in Fig. 1 include interactions 
of protons with 2.7 K CMBR. The synchrotron radiation of 
the secondary electrons peaks at GeV energies, with lumi- 
nosities at 1 keV of about 3.5 x 10^® erg/s for the case (a) 
and by a factor of 20 larger for the case (b) (not seen in 
Fig. 1). In the hot spot of Pictor A, the interactions of pro- 
tons with the local (synchrotron) radiation ai v ^ 0.1 — 100 
GHz are more (by two orders of magnitude) frequent, es- 
pecially if we assume a small, sub-kpc size of the hot spot. 
Even so, this process is not sufficiently effective to reproduce 
the observed X-ray fiuxes. Moreover, the synchrotron radia- 
tion of secondary (from jry interactions) electrons at 1 keV 
has, independent of the proton spectrum, standard spectrum 
with spectral index ~ 0.5, i.e. significantly harder than the 
spectrum detected by Chandra (Wilson et al. 2001). 

3.2 3C 120 

The X-radiation from a radio knot at a distance of 25" 
from the nucleus of the Seyfert galaxy 30 120 found in the 
ROSAT HRI data (Harris et al. 1999) reserves a special in- 
terest. The flux predicted by the SSC models appears several 
orders of magnitude below the observed X-ray flux. The two- 
component synchrotron model which postulates two differ- 
ent populations of relativistic electrons responsible for the 
radio-to-optical and X-ray emissions, at first glance seems 
a natural interpretation. However, the reported flux limit 
at optical frequencies (Harris et al. 1999) challenges this 
model as well. Indeed, the upper limit at f = 6.8 x 10^* Hz 
(Su < 0.18/iJy) combined with the X-ray flux at 2 keV 
(0.018/iJy) (see Fig. 1) requires an unusually hard spec- 
trum between the optical and X-ray band. Si, oc v'" with 
a < 0.35. This implies that the spectrum of radiating elec- 
trons should be flatter than N{Ec) oc E^^'"' . Although such 
a spectrum deviates from predictions of the canonical shock 
acceleration theory, formally we cannot exclude other, more 
effective particle acceleration scenarios. However, as long as 
the synchrotron energy losses [dEe/dt oc E^) play dominant 
role in formation of the electron spectrum, the latter cannot 
be harder, independent of the initial (acceleration) spectrum 
Q{Ee), than N{E^) oc E'^. Indeed, since the cooled steady- 
state spectrum of electrons is determined as 

N{E^) oc (dSe/di)-' / Q{E^)dEe, (14) 

even in the extreme case of abrupt low-energy cutoff in the 

injection spectrum (i.e. Q{Ec) = QqE~^° at Ec > E* , and 
Q{Ee) = at Ec < E"), the spectrum of cooled electrons 
has a broken power-law shape with N{Ee) oc E^^^"'^^^ at 
Ee> E*, and N{Ee) oc E'^ at E^ < E*. The acceleration 
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spectrum Q{Ea) can be sustained unchanged, assuming fast, 
with speed of light, escape or adiabatic losses of electrons. 
However, even with such a dramatic assumption, the elec- 
tron escape time, tcsc ~ R/c ~ 3 x 10^-Rkpc yr, in the knot 
with a radius 2-3 kpc appears longer than the typical lifetime 
of electrons responsible for synchrotron radiation of optical 
and X-ray photons, tsynch — 1.5-B~q ^(-B/l keV)"^''^ yr, un- 
less the magnetic field is less than ~ 10~^ G, or we should 
assume (Harris et al. 1999) that the observed X-ray flux is a 
result of contributions from short-lived (transient) compact 
regions where the energy-independent escape losses could 
dominate over the synchrotron losses. 

We do not face such a problem with the synchrotron ra- 
diation of protons. For a reasonable set of model-parameters, 
the proton escape could dominate over the synchrotron cool- 
ing, but yet the luminosity of the proton-synchrotron radi- 
ation could be maintained at a rather high level. Obviously, 
quite similar to the electron synchrotron model, the compar- 
ison of the X-ray and optical fiuxes tells us that the spectrum 
of the radiating protons should be harder than E~^''' . This 
implies that for even hardest possible acceleration spectra 
oc E~^'^ (Malkov 1999), the escape of particles should be es- 
sentially energy-independent (fcsc ex with P < 0.2). In 
Fig. 1 we show a possible version of the proton-synchrotron 
radiation calculated for the following model parameters: the 
magnetic field B — 1.5 mG, radius of the X-ray production 
region R — 3 kpc, acceleration spectrum of protons with 
Fo = 1.7 and Eo = 10^° eV, energy-independent escape of 
protons with tcsc = 5 x 10^ yr. For the diffusive propagation 
of particles, the latter corresponds to the diffusion coefficient 
which at proton energies > 10^^ (responsible for production 
of 0.1-10 keV synchrotron photons) is much larger than the 
Bohm diffusion coefficient. This provides the dominance of 
escape losses, and thus does not allow deformation of the 
primary (acceleration) spectrum of protons up to highest 
energies (see curve 3 in Fig. 2) . Then, the efficiency of the 
proton-synchrotron radiation, which is characterized by the 
ratio of the escape time to the synchrotron cooling time, at 
1 keV is about « 0.4 per cent. The efficiency is much higher 
(close to 100 per cent) for 100 MeV photons (see curve 3 in 
Fig. 1) produced by highest energy (~ 10'^'' eV) protons. The 
corresponding flux of 100 MeV 7-rays J^, ~ 5 x 10~* ph/cm^s 
is sufficiently large to be detected by GLAST (Gehrels & 
Michelson 1999). 

To reproduce the observed X-ray flux at 1 keV we must 
assume a rather high injection rate of accelerated protons 
during 10* yr operation of the jet, Lp = 1.65 x 10''^ erg/s. 
Note, however, that > 10^* eV protons do not contribute, for 
the given magnetic fleld of about 1.5 mG, to the production 
of < 1 keV photons. Therefore we may somewhat (by a 
factor of 4) soften this requirement by reducing the high 
energy cut-off in the proton spectrum down to Eq = 10^* 
eV. Further significant reduction of the required acceleration 
power could be achieved assuming more effective particle 
confinement. In Fig. 1 we show the luminosity of the proton- 
synchrotron radiation (curve 4) calculated for tcsc = 5 x 
10^ yiQ and Eo — 5x 10^* eV, but sustaining all other model 

^ In order to avoid violation of the Bohm limit at higher energies, 
it is assumed that above 10^* eV the escape time decreases as 
(X 1/E. 



parameters unchanged. In spite of the increase of the escape 
time by two orders of magnitude, which becomes comparable 
with the age of the source At = 10* yr, the synchrotron 
cooling remains a slower process, compared with the particle 
escape, until Ep — 1.2 x lO'^* eV. Therefore the break in the 
synchrotron spectrum occurs only at ~ 5 keV. At the same 
time, the required proton acceleration power becomes rather 
modest; now only Lp ~ 1.0 x 10*"^ erg/s is needed to explain 
the observed X-ray flux. 

3.3 PKS 0637-752 

This quasar at a redshift z = 0.651 has the largest and 
most powerful X-ray jet detected so far (Chartas et al. 2000, 
Schwartz et al. 2000). The X-ray luminosity of the > 100 
kpc jet of about Lx ~ 4 x 10** erg/s is contributed mostly 
by bright condensations between 7". 5 and 10". This region 
contains three 3 knots resolved in radio and optical wave- 
lengths. The enhanced X-ray emission associates, most prob- 
ably, with these knots, although its profile does not exactly 
repeat the radio and optical profiles of the jet. Therefore 
the combined flux of these three knots ~ 0.57 /iJy at 
an effective frequency 4.3 x 10^* Hz (Schwartz et al. 2000) 
should be considered as an upper limit for the optical flux 
from the X-ray emitting regions. 

Because the individual X-ray knots are not clearly re- 
solved, here we consider a simplified picture, namely treat 
the overall emission from the X-ray enhanced region as ra- 
diation from a single source. The X-ray luminosity based on 
the measured flux density at 1 keV, 5.9x 10~^* erg/cm^s keV 
and spectral index a = 0.85 (Schwartz et al. 2000) is shown, 
together with the model calculations, in Fig. 1. Curve 5 is ob- 
tained for the following parameters: B — 1.5 mG, Fo = 1.75, 
Eo = 10^° eV, i? = 5 kpc. Also it was assumed that par- 
ticles propagate in a "relaxed" Bohm diffusion regime with 
gyro-factor 77 = 10. Under these conditions the escape losses 
dominate over the synchrotron losses, and become important 
above 10^^ eV. This result in the steepening of the protons 
spectrum, from the initial _E~^'^^ to The latter con- 

tinues up to the exponential cutoff in the acceleration spec- 
trum at Eo = 10^" eV (curve 5 in Fig. 2). The corresponding 
proton-synchrotron spectrum nicely fits the observed X-ray 
spectrum and, at the same time, agrees with the optical flux 
reported from that part of the jet. Note, however, that the 
optical flux shown in Fig. 1 perhaps should be considered 
rather as an upper limit. Correspondingly, if the intrinsic 
optical flux of the X-ray emitting region is significantly less 
than this upper limit, we must assume flatter proton accel- 
eration spectrum. It should be noticed also that the extrap- 
olation of the proton synchrotron spectrum down to radio 
wavelengths appears significantly below the reported radio 
fluxes at 4.8 GHz and 8.6 GHz (Chartas et al. 2000), which 
most probably are due to the synchrotron radiation of di- 
rectly accelerated electrons. 

For the chosen set of model parameters, the absolute X- 
ray flux requires huge acceleration power Lp = 3x 10*® erg/s, 
which however agrees with the estimates of the jet power (in 
the form of kinetic energy of particles and of Poynting flux) 
which in most luminous extragalactic objects could be as 
large as 10*^ - 10** erg/s (see e.g. Sikora 2000, Ghisellini 
& Celotti 2001). We may nevertheless reduce the energy 
requirements assuming a somewhat different model param- 
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eters which would minimize the escape losses. An example 
of X-ray luminosity calculated for a favorable set of model 
parameters is presented by curve 6 in Fig. 1. In particular, 
compared with curve 5 the magnetic field is increased to 3 
mG, the size is reduced to 3 kpc, and it is assumed that 
the particle diffusion proceeds in the Bohm regime (r? = 1). 
This implies that B^gJi^pgr?"^ = 243, i.e. the proton losses 
are strongly dominated by synchrotron cooling. In order to 
reduce the required proton acceleration power by another 
factor of ~ 3, an "early" exponential cutoff in the accelera- 
tion spectrum is assumed at iSo = 2 x lO'^* eV. The resulting 
proton spectrum shown in Fig. 2 (curve 6) could be treated 
as an optimum spectrum designed to maximize the produc- 
tion of X-ray at 1 kcV, as it is seen in Fig. 1. The required 
proton power now is 10 times less than in the previous case, 
Lp = 2.9 X 10"*^ erg/s. Further significant reduction of Lp is 
almost impossible given the fact that the X-ray luminosity 
in the interval 0.2-10 kcV already is huge, Lx — 5 x 10*'' 
erg/s, unless wc assume that the jet is moving with a bulk 
Lorentz factor Fj ^ 1 towards the observer. 

Scarpa and Urry (1999) argued that the relativistic bulk 
motion of the jets on kpc scales, and the consequent beaming 
effect (if the jets are aligned with the line of sight) would 
allow reduction of the magnetic field, and consequently - 
longer electron lifetimes. This would help, at least partly, to 
remove some of the problems of electron synchrotron models 
suggested for radio and optical jets. Actually, the hypothesis 
of relativistic beaming with Doppler factor Dj ^ 1 implies 
dramatic reduction of the jet's intrinsic luminosity, and thus 
offers significantly relaxed conditions to all models of non- 
thermal emission of large scale extragalactic jets. In partic- 
ular, Tavecchio et al. (2000) and Celotti et al. (2001) have 
shown that the inverse Compton scattering of relatively low 
energy electrons on CMBR is a viable mechanism for expla- 
nation of X-ray emission of the jet of PKS 0637-752. Ob- 
viously, the beaming effect would be a great help also for 
the proton synchrotron model in the sense of significant re- 
duction of the required proton acceleration rate and/or the 
strength of the magnetic field. 



4 JET IN 3C 273 

At least 4 features have been revealed in the Chandra images 
of the X-ray jet of quasar 3C 273 (Marshall et al. 2000; Sam- 
bruna et al. 2001). Approximately 40% of the total X-ray lu- 
minosity of the jet is contributed by the brightest knot called 
Al (Bahcall et al. 1995). This knot is unique also in the sense 
that its radio, optical and X-ray fluxes are consistent with 
a single power-law spectrum (Rosner et al. 2000, Marshall 
et al. 2001; see, however, Sambruna et al. 2001). The spec- 
tral energy distribution (SED) of the knot Al is shown in 
Fig. 3. The observed fluxes at radio (Merlin array), optical 
(HST), and X-ray (Chandra) bands are from Marshall et all 
(2001). The overall slope with spectral index Qr-x — 0.75 
agrees well with the local slope at optical wavelengths, and 
matches the X-ray flux at 1 koV, although the local slope 
at X-rays is sliglitly different, Qx = 0.60 ± 0.05 (Marshall et 
all 2001). Rosner et al. (2000) has interpreted the power-law 
behavior of the broad-band spectrum of Al as a signature 
of synchrotron radiation by a single population of electrons. 
The spectral index 0.75 requires power-law index of elec- 




LogE (eV) 

Figure 3. The spectral energy distribution of the knot Al in 
3C 273. The radio, optical and X-ray fluxes are from Marshall 
ct al (2001). The solid, daslicd and dot-dashed curves corre- 
spond to 3 diff'crent sets of model parameters discussed in the 
text. Tlie heavy, standard and thin lines represent the fluxes of 
(1) synchrotron radiation of protons, (2) synchrotron radiation of 
secondary electrons produced in p7 and pp interactions, and (3) 
TF^-decay 7-rays produced at pp interactions, respectively. 

trons Fe = 2.5 and therefore, since the synchrotron cooling 
time of electrons is short compared with the lifetime of the 
jet, we should assume a hard, Q(Ee) oc E~^'^ type electron 
acceleration spectrum. Because of severe synchrotron losses, 
the X-ray emitting electrons could not propagate far from 
their birthplace/accelerator. To some extent this is true also 
for the electrons responsible for optical radiation. Therefore 
we should expect a point-source type morphology, unless the 
electron acceleration takes place throughout the entire knot, 
or the radiation is contributed by many compact regions in- 
side the knot. 



4.1 Proton-synchrotron model 

The proton synchrotron model offers more freedom for the 
principal model parameters. Indeed, for the size of the knot 
of about 1 kpc and for the ambient magnetic field B > 1 
mG, the synchrotron-cooling and the escape times of protons 
are comparable with the age of the jet. Thus, with certain 
assumptions about the acceleration spectrum and the prop- 
agation of protons we can satisfactorily fit the fluxes as well 
as to explain the kpc size of the observed diffuse nonthermal 
emission. 

In Fig. 3 we show three spectra of proton-synchrotron 
radiation modeled for the knot Al. For all 3 curves we as- 
sume that continuous injection of relativistic protons pro- 
ceeds with a constant rate during last 3 x 10^ yr. The solid 
curve is calculated for B = 5 mG and R = \ kpc, and assum- 
ing that the protons propagate in the Bohm diffusion regime 
with r? = 1 (model I). This implies that the synchrotron 
losses dominate over the particle escape, but yet, for the as- 
sumed source age and the magnetic field, the synchrotron 
losses become important only for protons responsible for X- 
ray flux above 1 keV {Ep > 3 x 10^*^ eV). Since the escape 
losses are also small and do not change the acceleration spec- 
trum of protons either, in order to explain the entire range 
of nonthermal emission from radio to X-rays with spectral 
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Figure 4. Total energy of relativistic protons trapped in the Icnot 
Al of the quasar 3C 273. The heavy solid, daslicd, and dot-dashed 
curves correspond to 3 different combinations of the model pa- 
rameters as in Fig. 3. The fluxes of secondary p7 and pp electrons 
multiplied by 10^" are also shown (thin lines). 



index Or-x — 0.7, we must assume a steep protou accelera- 
tion spectrum with To = 2Q;r-x + 1 = 2.4. The solid curve 
in Fig. 3 corresponds to such a power-law index of acceler- 
ated protons. The disadvantage of this fit is that it requires 
uncomfortably large acceleration rate ip = 1.2 x 10*^ erg/s 
with the total amount of protons deposited in the knot dur- 
ing 3 X 10'' yrs, ~ lO"^ erg (see Fig. 4). This imphes 
that the energy density of cosmic rays exceeds the density 
of the magnetic field by 3 orders of magnitude. 

It should be noticed, however, that this huge total en- 
ergy is contributed by low energy protons which, from the 
point of view of production of radiation above > 1 GHz, are 
in fact "vain" particles. Therefore assuming a low-energy 
cutoff in the acceleration spectrum, e.g. at 10^''' eV, we can 
significantly reduce the energy requirements. Another, per- 
haps a more natural way to avoid the low-energy protons can 
be achieved assuming a hard acceleration spectrum coupled 
with energy-dependent escape of particles. The dot-dashed 
curve in Fig. 4 corresponds to the case of the canonical shock 
acceleration spectrum with To = 2, and the escape time 
taken in the form te^dE) = 1.4 x 10^(£;/10" eV)"^/^ yr 
(model II). The resulting proton spectrum inside the knot 
in the most relevant energy region above 10^^ eV becomes 
steeper, E^^''' . The assumed large niaguctic field B = 10 
mG, as well as the upper energy cutoS' at i5 = 10^** eV 
in the acceleration spectrum allow additional reduction of 
the acceleration power down to 1.05 X lO'"^ erg/s. The cor- 
responding energy of protons confined in the knot now is 
4.9 X 10"" org. The spccti um of protou synchrotron radia- 
tion calculated for this set of parameters is shown in Fig. 3 
by the dot-dashed curve. 

Further significant reduction of the required en- 
ergy budget is still possible, provided that the proton- 
synchrotron radiation is responsible only for X-ray emis- 
sion, the radio and optical fluxes being related to other 
emission components, e.g. to the synchrotron radiation of 
electrons. The dashed curve in Fig. 3 corresponds to such 
a scenario (model 111). It is obtained assuming the follow- 
ing parameters: 5 = 3 mG, R = 2 kpc, Eq, = 10^* eV, 
Fo = 2, Lp = 10"'® erg/s. For these parameters the pro- 



ton synchrotron radiation explains the observed X-ray fiux. 
However, it appears below the measured radio and optical 
fluxes, because neither the synchrotron radiation nor the es- 
cape appear sufficient to cool effectively protons , and thus 
to steepen the original (acceleration) spectrum of protons 
at energies below < 10^^ eV (dashed curve in Fig. 4). A 
good compensation for such a compromise is the relatively 
modest (for 3C 273) injection power of relativistic protons. 

4.2 Synchrotron radiation of secondary electrons 

In Fig. 4 we show the spectra of secondary electrons pro- 
duced at interactions of protons with the ambient plasma 
and radiation flelds. The corresponding synchrotron radia- 
tion of these electrons is shown in Fig. 3. The energy den- 
sity of mm and sub-mm radiation, which is the most rel- 
evant band of electromagnetic radiation for interactions of 
highest energy protons, could be estimated from Fig. 3 as 
Wr = vSvidjH)^ jc K, 10~^^ erg/cm^, i.e. almost an order 
magnitude less than the energy density of CMBR at the 
epoch 2 = 0.158. 

The interactions of protons with CMBR do not con- 
tribute signiflcantly to the production of electrons with en- 
ergy less than 100 TeV, i.e. to the production of electrons 
responsible for synchrotron radiation at radio to X-ray en- 
ergies. The first generation electrons appear with energies 
exceeding Ec > lO^^eV. On the other hand, the large mag- 
netic field and small 7 — 7 optical depths prevent an ef- 
fective electromagnetic cascade development, which would 
allow production of lower energy electrons. As a result, the 
only signature of p7 interactions soon Fig. 3 is the fiatten- 
ing of the spectrum of the electron synchrotron radiation 
at McV/GcV energies (solid line), caused by electrons pro- 
duced through the Bethc-Hcitlcr pair production. 

For relatively high gas density in the knot, the interac- 
tions of accelerated protons with the thermal plasma con- 
tribute more efi^ectively, through decays of tt* mesons, in the 
production of relativistic electrons and positrons. Because of 
poorly known gas density in the jets, the calculations of sec- 
ondary "pp" electrons contain large uncertainties. On the 
other hand, an important information about the gas density 
could be provided by synchrotron radiation of these elec- 
trons; the latter obviously should not exceed the observed 
fluxes at radio, optical and X-ray bands. For the knots in 3C 
273, the most informative upper limits are contained in the 
very low radio fluxes at GHz frequencies. The curves in Fig. 3 
correspond to the gas density n = 10~® cm~^ for the model 
I, and 10~^ cm~^ for the models II and HI. It is seen that for 
the steep acceleration spectrum of protons oc E~^'* the radio 
data only marginally agree with the calculated synchrotron 
radio flux by secondary electrons (solid line in Fig. 3), and 
therefore n < 10~* cm~^. Such a strong constraint is a re- 
sult of very large amount of low, < 10 GeV energy protons 
(solid curve in Fig. 4) - the parents of secondary electrons 
producing synchrotron radio emission. Since the gas density 
in knots is likely to be signiflcantly larger, this constraint 
could be considered as an independent argument against 
the steep proton acceleration spectrum assumed in the the 
model 1. 

The models II and HI with flat E~^ type acceleration 
spectrum of protons allow much higher, respectively 10~^, 
and 10~^ cm~^ gas densities in the knot. Because these num- 
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Figure 5. Broad-band nonthermal emission of the knot Al pro- 
duced directly by accelerated protons via the synchrotron ra- 
diation (daslicd lines) and through the pp ^ tt'^ ^ f channel 
(dotted lines), and secondary electrons through the synchrotron 
(solid lines) and inverse Coinpton (dot-dashed lines) radiation. 
The heavy lines correspond to the case of slow escape (model 1) 
and thin lines correspond to fast escape (model II). The model 
parameters are discussed in the text. The synchrotron radiation 
of protons for the model II is too wealj and is beyond the figure 
frames. The contributions of TT^-decay 7-rays for the models I and 
II are almost identical, and therefore cannot be distinguished. 

bers are quite close to the densities expected in the environ- 
ments of large scale extragalactic jets, it seems an attractive 
idea of referring a fraction, or even the entire synchrotron 
spectrum of the knot to the synchrotron radiation of sec- 
ondary "pp" electrons. The production spectrum of these 
electrons with energies > 1 GeV coincides with the power- 
law spectrum of protons, thus it is quite easy to guess the 
power-law index of acceleration spectrum of parent protons, 
taking into account that for the magnetic field in the knot 
B ~1 mG and for the jet's age > 10® yr, the radio electrons 
effectively cool through the synchrotron radiation. Namely, 

To = 2 Qr-x ~ 1-5. 

The fluxes of synchrotron radiation of secondary "pp" 
electrons calculated in the coritirmous injection regime dur- 
ing last 10^ yr are shown in Fig. 5 for two combinations of 
model parameters: (a) i? = 1 kpc, B = 0.1 mG, ij = 1; 
and (b) R = 0.25 kpc, B = 0.03 mG, r? = 10. In both 
cases the same "power-law with exponential cutoff" accel- 
eration spectrum for protons is assumed with To = 1.5 and 
Eo = 10^* eV. These parameters imply essentially differ- 
ent escape times of protons with corresponding breaks at 
=i 3.5 X 10^® eV and 6.5 x lO" eV for the cases (a) and 
(b), respectively. The resulting breaks in the synchrotron 
spectra of secondary electrons appear at 1 MeV and 1 eV, 
respectively. It is seen the the case (a) explains quite well 
the observed nonthermal spectrum from radio to X-rays, 
while the case (b) can explain only radio and optical fluxes. 
Note that for chosen parameters the highest energy protons 
quickly escape the knot. This effect, combined with the low 
magnetic field, results in a dramatic drop of emissivity of 
the proton synchrotron radiation (dashed line in Fig. 5). 

The synchrotron radiation by secondary "pp" electrons 
in radiogalaxies and AGN is not a new idea. The basic 
problem of this hypothesis - the deficit of an adequate tar- 
get material - is clearly recognized since early 60's (see 



e.g. Burbidge, Burbidge and Sandage 1963). Our study of 
the Al knot in 3C 273 faces the same problem. In or- 
der to match the absolute fluxes of the observed nonther- 
mal radiation from knot Al, we must assume a very large 
product of the gas density and the proton injection power, 
nip = 8.5 X 10""^ erg/s cm . Thus for any reasonable accel- 
eration rate of protons Lp < lO*'^ erg/s, the plasma density 
in the knot should be close to 1 cm~^. Such high densities 
are not supported, however, by depolarization studies of the 
radio emission (see e.g. Roser et al. 2000). In principle, this 
problem could be essentially softened if the knot is moving 
towards the observer with a relativistic Doppler factor ^ 1. 



4.3 Relaxed model parameters due to relativistic 
beaming 

Recently, Sambruna et al. (2001) invoked the relativistic 
beaming effect to explain the X-ray emission from the knots 
in 3C 273 bj' the inverse Compton scattering of low en- 
ergy electrons, the main target photon field being the 2.7 
K CMBR. While this model can successfully explains the 
diffuse character of X-ray emission (the X-ray emitting elec- 
trons have long cooling times), it must assume a second, 
synchrotron component of radio to optical emission. Note 
that the electrons responsible for optical radiation should 
have very large energies of about 1 TeV, and correspond- 
ingly cool on very short times scales. Therefore they still 
require in situ acceleration (or re-acceleration) . 

Obviously, the hypothesis of relativistic beaming has 
a universally "positive" effect for all radiation models. In 
particular it helps to reduce significantly the energy re- 
quirements to both proton synchrotron and the secondary- 
electrons synchrotron models discussed above. Indeed, if the 
jet maintains the relativistic speed up to kpc scales, for the 
Doppler factor of about (5j ~ 5 (estimated for the inner 
jet of 3C 273; see e.g. Abraham & Romero 1999), these 
requirements become so relaxed that we can successfully 
explain the entire radiation of the knot Al merely by the 
proton synchrotron radiation. At the same time, the rela- 
tivistic beaming makes the interpretation of the nonthermal 
emission of the knot by secondary "pp" electrons an attrac- 
tive alternative with a quite reasonable product of the gas 
density, the jet lifetime , and the proton acceleration rate 
(n/0.1 cm-3) (AT/10'' yr) (Lp/lO-*^ erg/s) ~ 1. 

4.4 Gamma-ray emissivity of the knots 

The quasar 3C 273 is identified as powerful 7-ray emitter 
(Lichti et al. 1995) with a peak in the spectral energy dis- 
tribution at 1-10 MeV at the level of ~ 4 x 10""' erg/cm^s; 
the energy flux at 1 GeV is about ~ 2 x 10~^^ erg/cm^s. 
These fluxes exceed by several orders of magnitudes the flux 
of the proton synchrotron radiation (sec Fig. 3). Below we 
briefly discuss two other potential mechanisms of 7-radiation 
related to interactions of relativistic protons. 

Besides the cooling through synchrotron radiation, the 
secondary electrons release their energy also through inverse 
Compton radiation which may extend to very high energies. 
However, for the parameters favorable for the proton syn- 
chrotron radiation, in particular for the magnetic field B > 
1 mG, the contribution of the IC component is negligible. 



© 0000 RAS, MNRAS 000, 000-000 



Proton synchrotron radiation 11 



and therefore do not appear in Fig. 3. For smaller magnetic 
field, B < 0.1 mG, the IC fluxes are higher (Fig. 5), but still 
well below the sensitivity of 7-ray detectors. This is also 
true for "direct" 7r°-decay 7-rays produced at pp interac- 
tions. The fluxes of this component are shown in Fig. 3 and 
Fig. 5. Note that the sharp drops of 7-ray fluxes above 100 
GeV is caused by intergalactic photon-photon absorption. 
In calculations presented in Fig. 3 and 5, we adopted one 
of the recent models of the difi^use extragalactic background 
of Primack ct al. (2001) which satisfactorily describes the 
observational data at near infrared and optical wavelengths 
- the most important band from the point of view of inter- 
galactic absorption of 7-rays above 100 GcV. 

Although the flux of 7r"-dccay 7-rays is proportional 
to the ambient gas density n and the total amount of ac- 
cumulated relativistic protons, Wp ~ AtLp, it cannot be 
arbitrarily increased assuming larger values for the prod- 
uct riAtLp. The 7r°-decay 7-rays arc tightly coupled with 
the synchrotron radiation of secondary (vr^-decay) electrons. 
Therefore the fluxes of 7r°-decay 7-rays are (unavoidably) 
limited by the energy fluxes of synchrotron radiation of sec- 
ondary electrons at low frequencies. In particular, since both 
GcV 7-rays and GHz radio photons arc initiated by the same 
(10 to 100 GeV) primary protons, the flux of 7-rays around 
1 GeV cannot exceed 10~^^ erg/cm^s, as it is clearly seen 
in Fig. 3 and Fig. 5. Obviously this upper limit is insen- 
sitive to the relativistic beaming cfi^ects. This excludes any 
chance to explain the observed high energy 7-radiation from 
3C 273 by pp interactions in the the large-scale jet. The ob- 
served MeV/GeV 7-radiation originates, most probably, in 
the inner jet. 



5 RUNAWAY PROTONS OUTSIDE JETS 

Beyond 100 Mpc, the Universe is opaque for extremely high 
energy (EHE) cosmic rays with energy E > 10^° eV (see 
e.g. Cronin 1999). Fortunately, observations of the electro- 
magnetic radiation associated with these particles can com- 
pensate, at least partly, this limitation, allowing extension 
of the (indirect) study of EHE cosmic rays to cosmological 
scales. 

The results of previous sections show that an essen- 
tial fraction of the nonthermal energy of the jet releeised in 
the knots and hot spots in the form of relativistic protons, 
eventually escapes the jet; the kpc scales of these structures 
are not sufficient for effective confinement of particles with 
E > 10^® cV. Thus the nonthermal radiation of the jet docs 
not tell us much about the most energetic particles acceler- 
ated in the knots and hot spots. Such a crucial information 
can be recovered by observing characteristic nonthermal ra- 
diation of the cluster environments harboring the runaway 
particles. The spatial and spectral characteristics of non- 
thermal emission initiated by interactions of protons with 
the CMBR and the ambient gas, generally depend on both 
the spectrum of runaway protons and the properties of the 
local environment. Below I will discuss two different cases, 
assuming that (1) the source of EHE protons is surrounded 
by a rich galaxy cluster with B-field exceeding 1 /xG, and 
(2) the EHE source is located in a low magnetic field envi- 
ronment. 



5.1 Radiation initiated in rich galaxy clusters 

It is believed that a non-ncgligiblc fraction of the energy 
of clusters of galaxies is contained in non-thermal forms, 
in particular in the form of relativistic particles and mag- 
netic fields. While this conclusion concerning the energy 
budget of relativistic particles is still based on theoretical 
and phenomcnological arguments comprehensively discussed 
30-years ago by Brecher & Burbidge (1972), there is grow- 
ing evidence that many rich galaxy clusters have magnetic 
field strengths at the microgauss or even higher level (Kro- 
nbcrg 2001). Remarkably, this concerns not only the rich 
galax;y clusters with strong cooling-flows, but also the "nor- 
mal", non-cooling- fiow clusters (Clarke et al. 2001). The es- 
timated magnetic field of about 5 fxG with energy density 
wb = B^/8tv « 0.6 eV/cm^ in the inner _R ~ 0.5 Mpc 
sphere ( Clarke et al. 2001), would result, most probably, 
in relatively slow diffusion of cosmic rays in the intracluster 
medium, and thus would provide their effective confinement 
and accumulation in the galaxy clusters (Volk et al. 1996, 
Berezinsky et al. 1997, Blasi & Colafraricesco 1999). 

The cnergy-dcpcudcut diffusion coefficient, together 
with the initial spectrum of particles injected in the intr- 
acluster medium, determine the spectrum and the total en- 
ergy content of protons established in the cluster over life- 
times of about t ~ 1/Ho ~ 10^" yr. In Figs. 6a and 6b 
we show two examples of contemporary proton spectra ob- 
tained within the "leaJjy-box" type approximation of parti- 
cle propagation (see e.g. Ginzburg & Syrovatskii 1964) for 
the following model parameters: 

(a) continuous injection of relativistic protons with a con- 
stant rate Lp — 10^ ' crg/s during At = 10^° yr; acceleration 
spectrum of protons with power-law index Fp = 2 and ex- 
ponential cutoff at Eq = 10^° eV; intracluster magnetic field 
B = 3 nG and gas density n = 10~^ cm~^ within the central 
R — 0.5 Mpc region; escape time of protons from the high 
magnetic field region Tcsc = 5 x lO^iJj,/'''^ yr; 

(b) Lp = 3 X 10""* erg/s; At = 10* yr; Fp = 1.5, Eq = 
10^" cV, R = 0.5 Mpc, B = 5 /iG, n = 10-=^ cm'^, Tesc = 
2.0 X lO^Sj^/''^ yr. 

For the assumed magnetic field B and the cluster radius 
R, the chosen escape times correspond to diffusion coeffi- 
cients given in the form D{E) = DoE^g with Do and /3 com- 
patible with the model parameters discussed in the litera- 
ture (sec e.g. Volk ct al. 1996, Berezinsky et al. 1997, Blasi & 
Colafrancesco 1999). For the given diffusion coefficients, the 
particle escape becomes an important factor in formation of 
the proton spectra above E* oc (R^ /Do At^^^ 3x lO" eV 
and lO^'^ eV, for the cases (a) and (b), respectively. As a 
result, respectively E~'^'^ and Ep^'^^ type proton spectra 
are formed above these energies, before approaching to the 
intrinsic (acceleration) cutoff at E ^ 10^" eV. 

The cases (a) and (b) correspond to two essentially dif- 
ferent scenarios. The case (a) could be treated as a quasi- 
continuous injection of protons into the cluster during its 
age of about 10^° yr. In this case not only AGN jets, but 
also other sources of particle acceleration, e.g. large-scale 
(multi-hundred kpc) shock structures (Miniati et al. 2001) 
can contribute to the high energy proton production. The 
case (b) correspond to the scenario when the proton produc- 
tion is contributed by a single powerful AGN (see e.g. En- 
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sslin et al. 1997) over its lifetime of about 10* yr. It should 
be noticed that, because the escape time of protons with 
energy E > 10^* eV is less than 10* yr, the relatively short 
time of operation of the accelerator has an impact only on 
the accumulated energy budget of lower energy protons. In 
other words, the amount of highest energy protons confined 
in the cluster is effectively determined by recent accelera- 
tor(s) operating since last 10^ — 10* years. Also, the weak 
energy dependence of the escape time (oc E~^^^) and the 

very hard proton spectrum (oc with the injection 

rate, Lp = 3 x 10"*^ erg/s, provide much higher content of 
> 10^* eV protons confined in the cluster for the case (b) 
compared with the case (a) (see Fig. 6). The assumed pro- 
ton acceleration rate is large, but still acceptable, given the 
fact that the total kinetic energy of electron-proton jets in 
powerful radio sources can significantly exceed 10'*'^ erg/s 
(Ghisellini & Celotti 2001), provided that the kinetic energy 
power at the core is conserved along the jet up to > 100 kpc 
scales (Celotti & Fabian 1993). The numerical calculations 
show that, for the given duration and rate of proton injec- 
tion, the total energy contents of protons accumulated in 
the cluster are Wp = 1.4 x 10^^ erg and 2.1 x 10^^ erg for 
the cases (a) and (b), respectively. These numbers do agree 
with general phenomenological considerations discussed in 
the literature (Brecher and Burbidge 1972, Volk et al. 1996, 
Berezinsky et al. 1997, Ensslin et al. 1997, Lieu et al. 1999). 

In Fig. 6 we show the spectra of secondary electrons 
produced at interactions of relativistic protons with the in- 
tracluster gas for n = 10"^ cm"^, and with 2.7 K CMBR. 
These spectra consist of three components: (i) electrons from 
pp interactions, (ii) electrons from p'y interactions due to 
the Bethe-Heitler pair production and (iii) electrons from 
the photo-meson processes. Due to the radiative cooling, 
which in the case of electrons strongly dominates over es- 
cape losses at all energies, the first population electrons have 
a power-law spectrum oc Ec while the spectra of the 

second and third electron populations are almost indepen- 
dent of the spectrum of primary protons, namely both cooled 
electron populations have standard E~^ type spectra with 
high-energy cutoffs at energies ~ 10^^ eV and ~ lO'^^ eV, re- 
spectively. While in the case (a) the second and third (p7) 
electron components are strongly suppressed (Fig. 6a), in 
the case (b) the features of all three components are clearly 
seen in Fig. 6b. 

Although the energy content in secondary electrons is 
very small compared with the total energy budget, approxi- 
mately half of the nonthermal energy of the cluster is even- 
tually radiated through these secondary electrons. It is re- 
markable, that unlike primary (directly accelerated) elec- 
trons, which because of short lifetimes are concentrated in 
the proximity of their accelerators, the secondary electrons 
are homogeneously distributed over the entire cluster, and 
therefore their radiation has an extended (diffuse) character. 

The production rates of nonthermal radiation consist- 



I Note that due to the confinement of low-energy particles in the 
jet, the spectrum of protons injected in the intracluster medium 
should have a low-energy cutoff. But since this effect does not 
have a noticeable impact on the results, and in order to avoid 
unnecessary increase of model parameters, we assume a single 
power-law injection spectrum. 
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Figure 6. The total energy content of protons and secondary 
electrons (multiplied by 10^) inside a galaxy cluster calculated 
for scenarios (a) and (b) (see Sec. 5.1 for the assumed model 
parameters). 



ing of 5 components - (1) synchrotron (marked as e-synch) 
and (2) inverse Compton (IC) photons emitted by secondary 
electrons, (3) synchrotron radiation of protons (p-synch), 
and 7r°-decay 7-rays from (4) proton-proton (pp) and (5) 
proton-photon (p-y) interactions - are shown in Figs. 7a and 
7b. It is seen that the spectral energy distributions (SED) 
of radiation characterizing the (a) and (b) scenarios are 
essentially different. 

In the case (a) the nonthermal radiation is mainly con- 
tributed, directly or via secondary electrons, by pp interac- 
tions . The importance of these interactions in the clusters 
of galaxies has been often discussed in the literature (Ves- 
trand 1982; Dermer & Rephaeh 1988; Dar & Shaviv 1995; 
Volk et al. 1996; Berezinsky et al. 1997; Ensslin et al. 1997; 
Atoyan & Volk 2000). Since the pp interaction timescales in 
the galaxy clusters with n < 10"'^ cm~'^ exceed the source 
ages of about 10^° yr, the absolute fluxes of radiation are 
proportional to the product nLpAt. For the assumed in- 
dex of accelerated protons Fq — 2, approximately the same 
fraction of the proton kinetic energy is released in yr'^-decay 
7-rays and TT^-decay electrons and positrons. On the other 
hand, the assumed magnetic field B = 3 fiG implies an en- 
ergy density close to the density of 2.7 K CMBR, therefore 
equal fractions of the electron energy are released through 
the synchrotron and inverse Compton channels. This results 
in the flat overall SED over a very broad frequency range 
from radio to multi-TeV 7-rays (Fig. 7a). The energy domain 
below 1 keV is due to synchrotron radiation of electrons, 
while the interval between X-rays and low energy (< 100 
MeV) 7-rays is contributed by inverse Compton mechanism. 
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Figure 7. Broad-band luminosities of nonthermal radiation ini- 
tiated by protons in a galaxy cluster. The calculations correspond 
to two different scenarios of proton injection into the cluster (see 
the text for details). 
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Figure 8. Broad-band SED of nonthermal radiation of clusters 
corresponding to luminosities shown in Fig. 7a and Fig. 7b, and 
assuming that the source is located at d = 175 Mpc (a) and 790 
Mpc (b), respectively (see the text for details and used model 
parameters. 



At higher energies the radiation is dominated by 7r°-decay 
7-rays produced at pp interactions. The local maximum at 
10^^ — 10^" eV in Fig. 7a, due to decays of vr^-photomesons, 
is relatively weak because of effective escape of protons from 
the cluster (see Fig. 6a) 

For harder spectra of accelerated protons, e.g. with 
power-law index Fp — 1.5 the SED is strongly dominated 
by synchrotron radiation of secondary p7 electrons with two 
prominent peaks at MeV and TeV energies^ (Fig. 7b) cor- 
responding to the radiation by electrons from Bethe-Heitler 
pair production and photo-meson production processes, re- 
spectively. While at optical and radio frequencies the energy 
flux decreases as {uSi, oc i''^'^), at extremely high energies 
E 10^® — 10^° eV the prominent "p7" peak dominates over 
the entire SED (Fig. 7b). However, this peak, as a part of the 
whole region of 7-rays above 10 TeV is not visible for the ob- 
server. Because of interactions with the diffuse extragalactic 
background radiation, these energetic 7-rays disappear dur- 
ing their passage from the source to the observer. 



** Note that the synchrotron radiation by primary (directly ac- 
celerated) electrons has an unavoidable, self-regulated cutoff be- 



low |a~^mec2r;-l ~ 160 MeV (Aharonian 2000), the latter be- 
ing determined by the balance between the synchrotron loss rate 
and the maximum-possible (t; = 1) acceleration rate. Obviously 
in the case of synchrotron radiation of secondary electrons, there 
is no intrinsic limit on 7-ray energy. 



5.1.1 Detectability of gamma-rays 

In Fig. 8 we show the expected fluxes of radiation assum- 
ing that the sources with luminosities presented in Figs. 
7a and 7b are located at small (like Pictor A) 175 Mpc 
(case a) and large (like 3C 273) 790 Mpc (case b) distances. 
These are two specific examples chosen to have a feeling for 
the detectability of radiation at different wavelengths. For 
the given source luminosities, the fluxes of high energy 7- 
rays are determined not only by the distance to the source 
{Fry (X but also by the intergalactic photon-photon 

absorption. Note that the mean free path of 7-rays above 
100 TeV is less than 1 Mpc, so strictly speaking we should 
include in calculations the radiation of next generation elec- 
trons produced inside the cluster. However, since the 7-ray 
luminosities sharply drop above 100 TeV (see Fig. 7), we can 
ignore these (3rd generation) photons without a significant 
impact on the accuracy of calculations. Such an approxima- 
tion, however, could be inappropriate for 7-rays with energy 

> 10^^ eV. The mean free paths of these 7-rays, interacting 
with CMBR and extragalactic radio background photons, 
contain large uncertainties because of the lack of reliable 
information about the radio background at MHz frequen- 
cies (Berezinsky 1970, Aharonian et al. 1992, Protheroe & 
Biermann 1996, Coppi & Aharonian 1997). Nevertheless it is 
likely that the extragalactic photon fields cannot prevent the 

> 10'^ eV 7-rays to travel freely over distances more than 1 
Mpc. Actually the interaction of the highest energy 7-rays 
with the cluster's own synchrotron radio emission could be 
more important process. It is easy to show that the lO'^® eV 
7-rays can escape the central 1 Mpc region of the clus- 
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ter if the radio luminosity at frequencies ~ 10 MHz docs 
not exceed 10*^erg/s. For comparison, the 10 MHz radio 
luminosity of the Coma cluster is about 10*" erg/s, but in 
most powerful clusters it could be significantly higher. For 
example, the power-law extrapolation (with spectral index 
1.6) of the recently reported radio flux from the giant halo 
in Abell 2163 (2 = 0.203) at 20 cm ~ 0.15 Jy, to lower 
frequencies (Feretti et al. 2001) gives for the luminosity of 
the cluster ~ 5 x 10*^ erg/s at 10 MHz. This impHes that 
the > lO^^oV 7-rays arc absorbed inside this cluster. In 
such cases we have to include in the overall SED the syn- 
chrotron radiation of the secondary electrons. This would 
simply result in the linear increase (by a factor of two) of 
the flux of > 1 TeV 7-rays (typical energy of synchrotron 
photons produced by W^^ eV electrons in the magnetic field 
B > 1 nG) shown in Fig. 8b. Even so, it would be difficult 
to detect this radiation due to the heavy intergalactic ab- 
sorption of > 0.1 TcV 7-rays arriving from cosmologically 
distant (d 3> 100 Mpc) sources (see e.g. Steckcr ct al. 1992). 
The chances to detect the "echo" of the original > 10^® eV 7- 
rays would be dramatically increased if the free path of these 
photons exceeds 1 Mpc. In this case the photon-photon inter- 
actions take place predominantly outside the cluster where 
the intergalactic magnetic field Big is significantly smaller 
(most probably. Big < 10~* G). Then we should expect an 
extended synchrotron emission from the pair-produced (and 
radiatively cooled) electrons with hard differential spectrum 
oc E~^'^ up to the maximum energy in the SED expected at 

Sn,ax ~ 5(Big/10-^ G) (Ee/lO'^ eV)^ GeV . (15) 

The 7-ray spectrum beyond -Emax depends on the shape of 
the proton spectrum in the region of the cutoff Eo, but 
in any case it is smoother (a basic feature of the syn- 
chrotron radiation) than the proton spectrum in the cor- 
responding region above Eq. The energy flux of radiation 
which peaks at iJmax can be easily estimated from Fig. 7b, 
assuming that almost the entire 10'*'' erg/s luminosity in 
the primary E ~ 10^® — 10"^" cV 7-rays is re-radiated out- 
side the cluster in the form of secondary synchrotron pho- 
tons. Namely, we should expect an energy flux at the level 
10^"*^^ (d/1 Gpc)^'^ crg/cm^s contributed mainly by 
the region around -Bmax. It is interesting to note that al- 
though the secondary synchrotron radiation is produced out- 
side the cluster, the angular distribution of this radiation 
seen by the observer would essentially coincide with the size 
of the cluster. Indeed, the ratio of the gyroradius to the 
mean synchrotron interaction path of electrons, which is es- 
timated as r-g/Asynch — 30 EfglBic / 10~^ G), implies that 
the photo-produced electrons of energy Ec ~ 10^^ eV will 
radiate synchrotron photons before they would significantly 
change their direction in the mag netic field > 10"® G, i.e. the 
synchrotron 7-rays produced outside the cluster will follow 
the direction of their "grandparents" - 7r°-decay 7-rays pro- 
duced inside the cluster. Below -Emax the angular distribu- 
tion of 7-rays could be quite different. In this energy region 
7-rays produced by cooled electrons with steady state dis- 
tribution oc E~'^, have hard power-law differential spectrum 
oc S"^-^. For Big ~ 10"® G, the trajectories oi E < 10^* eV 
electrons will be curved significantly before they radiate syn- 
chrotron 7-rays. Therefore, the size of the low energy 7-rays, 
E < 100 MeV, produced by these electrons, is determined 
by the mean free path of primary ~ 10^® eV 7-rays in the 



intergalactic photon field, and thus would exceed the an- 
gular size of the cluster. The performance of the GLAST 
instrument (see e.g. Gehrels & Michelson 1999) is perfectly 
suited to both the expected fluxes and the energy range of 
this radiation even from cosmologically distant (rf > 1 Gpc) 
sources. Detection of such radiation would not only conflrm 
the presence of highest energy protons in the clusters of 
galaxies, but also could provide unique information about 
the intergalactic magnetic fields and the diffuse extragalac- 
tic radio background. 

5.1.2 Detectability of X-rays 

The "hadronic" nonthermal radiation luminosities shown in 
Fig. 8 depend on several model parameters like the mag- 
netic field and the gas density of the intracluster medium, 
the total amount of protons confined in the cluster, etc. The 
shape of radiation strongly depends on the spectrum of in- 
jected protons, as well as on the diffusion coefficient. How- 
ever, since both the assumed ambient gas density and the 
proton injection rates are close to the maximum allowed 
numbers, the luminosities shown in Fig. 8 hardly could be 
dramatically increased. If so, the nonthermal X-ray fluxes of 
"hadronic' origin could be marginally detected only from 
relatively nearby clusters like Coma, Perseus, the cluster 
surrounding Cygnus A, etc. The case of the Coma cluster 
presents a special interest because of hard, most probably 
nonthermal X-radiation recently reported from this source 
(see e.g. Rephaeli 2001). The most natural interpretation of 
this radiation by the inverse Compton scattering of radio- 
electrons requires an intracluster magnetic field ~ 0.1 /iG 
(e.g. Atoyan & Volk 2000), i.e. noticeably smaller than the 
the value of about several /iG deduced from Faraday ro- 
tation measurements. Although this discrepancy perhaps 
should not be overemphasized (for possible explanations see 
Rephaeli 2001, and Petrosian 2001), it is worth exploring 
other possibilities, in particular invoking the radiation of 
secondary electrons produced at interactions of EHE pro- 
tons with the ambient gas or photon fields. 

Concerning the scenario (a). Fig. 7a predicts a radio 
luminosity ~ 10 erg/s which is by two orders of magni- 
tude higher than the luminosity observed from Coma. Be- 
cause of the assumed large magnetic field B = 3fj,G, we 
face the same problem which arises in the standard inverse 
Compton models. This imposes an upper limit on the X- 
ray luminosity, Lx < 10''° erg/s. However, this should not 
be treated as a robust constraint. In the "hadronic" models 
this limit can be easily removed by assuming harder (than 
E~^) power-law proton spectrum, and/or a low-energy spec- 
tral cutoff. Of course, this would automatically suppress also 
the inverts Compton component of radiation produced by 
low energy secondary electrons. In contrast to the conven- 
tional inverse Compton models, in the "hadronic" model 
we have an additional important channel for production of 
X-rays - synchrotron radiation of very high energy (multi- 
TeV) secondary electrons. This is demonstrated in Fig. 7b 
(scenario (b)). Even so, for the distance to Coma ~ 100 
Mpc, the expected hard X-ray flux is below, by an order of 
magnitude, the reported 20-80 keV flux from Coma of about 
2x 10"^^ erg/cm^s. In the "hadronic" model the synchrotron 
X-ray luminosity slightly depends on the strength of the 
magnetic field, thus in order to compensate this deficit, we 
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have to increase the proton injection rate by the same fac- 
tor, i.e. up to > 10^^ erg/s. This uncomfortably high, at least 
for an object like Coma, acceleration power makes problem- 
atic, although cannot completely exclude, the "hadronic" 
origin of X-rays. Decisive tests could be provided by new 
spectroscopic measurements (the hadronic model unambigu- 
ously predicts a hard X-ray spectrum with photon index 
~ 1.5 up to 1 MeV). Also, the model predicts large 7-ray 
fluxes both at MeV/GeV and TeV energies. These compo- 
nents have, however, different origins. The radiation up to 
10 GeV is produced by "Bethe-Heitler" electrons, while the 
TeV 7-rays are due to the electrons from tt* decays. While 
the TeV luminosity would be suppressed if the exponential 
cutoff in the proton spectrum occurs below 10^° eV, with- 
out any impact on the X-ray flux, the MeV/GeV 7-rays are 
tightly connected with X-rays. This implies that we should 
expect similar fluxes in < 100 keV X-rays and > 100 MeV 
7-rays, i.e. if the "hadronic" origin of X-rays is correct, then 
F^{> 100 MeV) ~ 10"" erg/cm^s. This hypothesis can be 
easily checked by GLAST. 

Finally we note that the detection of "hadronic" X-rays 
from extended regions of clusters of galaxies is a rather hard 
task, even for instruments like Chandra and XMM-Newton, 
especially because of the presence of high local X-ray com- 
ponents like the thermal X-ray emission of the hot intra- 
cluster gas, nonthermal X-rays due to inverse Compton ra- 
diation of directly accelerated electrons, etc. Regardless of 
the details, the spectral band of high energy 7-rays seems a 
more promising window to explore the "hadronic" processes 
with GLAST and, perhaps also with forthcoming 100 GeV 
threshold Cherenkov telescope arrays. 

5.2 Very high energy radiation produced in 
low-magnetic-field environments 

The characteristics of nonthermal radiation of runaway pro- 
tons in the case of a lack of strong cluster environment 
surrounding the EHE source, are essentially different from 
the radiation features described in the previous section. 
The weak magnetic fields in such regions make faster the 
particle propagation, as well as prevent the dramatic syn- 
chrotron cooling of highest energy (first generation) elec- 
trons. This allows an effective development of relativistic 
electron-photon cascades triggered by interactions of run- 
away protons with 2.7 CMBR. The first stage of the cascade 
initiated by secondary electrons and 7-rays interacting with 
the same 2.7 K CMBR leads to formation of a standard 7-ray 
spectrum which can be approximated as dN^/dE oc E~^'^ 
at E < 10 TeV, and dN^/dE oc S'^-^^ at E > 10 TeV 
with a sharp cutoff at £ ~ 100 TeV. After the fast devel- 
opment in CMBR, the cascade enters the second (slower) 
stage. At this stage E < 100 TeV 7-rays produce pairs on 
the IR/O diffuse background radiation, while the Compton 
scattering of electrons is still dominated by CMBR photons 
(Berezinsky et al. 1990; Protheroe & Stanev 1993; Coppi 
& Aharonian 1997). The second-stage cascade, which actu- 
ally consists of 2-3 interactions, shifts the spectrum to lower 
(TeV and sub- TeV) energies, and broadens the angular dis- 
tribution of emission due to the deflections of iJ < 100 TeV 
electrons in the ambient magnetic field. Unfortunately, the 
intergalactic magnetic fields and their fluctuations in very 
large (^ 1 Mpc) scales remain highly unknown, which does 



not allow us to make definite conclusions concerning the ex- 
pected characteristics of radiation. Nevertheless, depending 
on the strength of intergalactic magnetic field, one of the 
following 7-ray emission components can be predicted. 

Very week magnetic field. For a very weak intergalac- 
tic field, B < 10^^^ Gp^ the cascade radiation arrives, 
because of almost rectilinear propagation of primary pro- 
tons and secondary pairs, from a direction centered on 
the source. In this case we may expect a point-like source 
of radiation, although the 7-rays are produced at dis- 
tances > 10 Mpc from the source. Indeed, for the given 
energy of a detected 7-ray photon EtoV = E/1 TeV, 
the emission angle is determined by the direction of elec- 
trons participating in the last interaction, namely by the 
deflection in the magnetic field of the parent electron 
of energy Ec = {E / AkTf/^m^c^ ~ 17£tcV TeV. The 
mean attenuation path of these electrons in the CMBR 
is about Ae — G.Q2Erj,^y Mpc, while the gyroradius is 
rg ~ 20 St{v(-Big/10"^'^G)-^Mpc. Correspondingly, the 
emission angle 61(e) ~ A^/rg ~ 10'^£.^ev (-Big/10"^^ G). 
Thus, for the intergalactic magnetic fleld of about 10~^® G, 
the cascade radiation at E > 100 GeV would be concen- 
trated within an angle of 1°. Since for sources at distances 
between 100 Mpc and 1000 Mpc, we expect hard cascade 
spectrum with a cutoff between 100 GeV and 1 TeV, an ap- 
proximately half of the energy of > 10^° eV protons (com- 
pletely lost at interaction with CMBR over distances < 100 
Mpc) would be released in this energy interval, and give a 
flux 

~ 5X10- f^^i^) erg/cm^s(16) 
^ \^ 10« erg/s y \^ 1 Gpc y ^' ^ ' 

Down to lower energies, the energy flux decreases (oc 
E^^'^), and the angular size of emission increases (oc E). 
Therefore searches for such an emission from directions of 
nonthermal extragalactic sources, in particular from AGN 
with powerful X-ray jets, can be done most effectively by 100 
GeV threshold imaging atmospheric Cherenkov telescopes 
(e.g. Aharonian & Akerlof 1997). 

Intermediate magnetic field. If Big > 10~^^ G, the cas- 
cade electrons are promptly isotropized. This leads to the 
formation of giant pair "halos" surrounding strong extra- 
galactic TeV sources (Aharonian et al. 1994). The angular 
size of the extended 7-ray source depends on photon en- 
ergy. For the given energy of the detected photon -BtgV , it is 
mainly determined by the mean free path of previous gen- 
eration 7-rays of energy E' ~ 2Ec ~ SAE^^y TeV. The free 
path of _B' ~ 10 — 15 TeV photons, which are responsible 
for the detected (last generation) 100 GeV cascade 7-rays, 
presently is poorly known, but, probably, it does not exceed 
50 Mpc (see e.g. Primack et al. 2001). Thus the typical size 
of a 100 GeV halo radiation surrounding the EHE source 
at a distance 1 Gpc would be < 3°. The detection of pair 
halos presents a difficult experimental task, compared, in 

tt Although quite speculative, such a large-scale intergalactic 
magnetic field cannot be a priori ruled out (see e.g. Plaga 1995; 
Waxman & Coppi 1996), in particular if an essential fraction of 
the Universe consists of huge, 100 Mpc scale voids (Einasto 2001). 



© 0000 RAS, MNRAS 000, 000-000 



16 F.A. Aharonian 



particular, with the detection of rectilinear cascade radia- 
tion discussed above. At the same time, as a compensation, 
the EHE sources can be revealed through their halo radiar 
tion independent of the orientation of AGN jets. It should 
be noticed in this regard that only in the case of an isotropi- 
cally emitting source the halo will be centered on the source. 
If the relativistic outflow injecting EHE protons is directed 
away from the observer, the center of the halo would be dis- 
placed by an angle comparable to the typical angular size 
of the halo. The radiation characteristics of halos initiated 
by EHE protons primarily depend on the level of the dif- 
fuse extragalactic background, first of all at mid- and far- 
infrared wavelengths, but not on the intcrgalactic magnetic 
field, provided that the latter does not exceed 10^^ G. 



Large magnetic field. When considering electromagnetic 
cascades initiated by EHE protons in the intcrgalactic 
medium with magnetic field Big > 10~® G, we must distin- 
guish between two populations of secondary electrons. The 
electrons originating from the Bethe-Heitler pair production 
process are produced with typical energies {me/mp)Ep ~ 
IQis _ iQis eV. For magnetic field < 10"^ G they cool 
mainly through inverse Compton scattering, and thus pro- 
duce faint (e"'",e~) halos in a way discussed above. Mean- 
while, the electrons originating from the photo-meson pro- 
duction process, directly, via tt "'"-decays, or through interac- 
tions of 7r°-decay 7-rays with the extragalactic difi'use radio 
background, have much higher energies, ~ l/lOSp > 10^^ 
eV, taking into account that only > 10^° eV protons in- 
teract effectively with CMBR. Because of the Klein-Nishina 
effect, the energy losses of those electrons are dominated 
by the synchrotron radiation, as long as the ambient mag- 
netic field exceeds ~ 10"^ G (Gould & Rephaeh 1978). 
This prevents the cascade development, but instead pro- 
vides another effective channel for production of high en- 
ergy 7-rays. Almost the whole energy of > 10^" eV pro- 
tons is released, through the synchrotron radiation of sec- 
ondary electrons, into the 7-rays with characteristic energy 
tmax ~ 50(5/10"** G)(E(,/10^" eV)^ GeV. Because the gyro- 
radius of lO'^" eV protons in the magnetic field Big > lO^'' G 
is comparable or less than their mean p7 interaction path 
Ap ~ 100 Mpc, we should expect a diffuse radiation compo- 
nent emitted by huge intcrgalactic regions. This diminishes 
the chances for detection of this radiation component, es- 
pecially from relatively nearby (d ^ 1 Gpc) objects. The 
situation is quite different for cosmologically distant ob- 
jects. Because at distant cosmological epochs the CMBR 
was denser, riph oc (1 -I- z)^, and hotter, Tr oc I + z, the 
mean free path of protons Ap has a strong z-depcndence. 
At energies _Bp < 3 x 10^" eV it can be approximated as 
Ap ~ 5.2(1 + zfex.p[3 x 10^° eV/(l + z)E] Mpc (Berezin- 
sky & Grigoreva 1988). For example, in the environments of 
quasars 3C 273 (z = 0.158) and PKS 0637-752 (z = 0.651), 
the mean free paths of lO^'' eV protons are 44.6 Mpc and 7.1 
Mpc, respectively. This implies that, if the X-ray emission 
from large scale jets of these powerful quasars has indeed 
proton-synchrotron origin, we may expect an accompanying 
GeV 7-radiation component initiated by highest energy run- 
away protons outside the jets, but still within ~ 3° centered 
on 3C 273, and within 10 arcminutes centered on PKS 0637- 
752. For some model parameters, discussed in Sections 3 and 



4, the expected fluxes of this component of radiation would 
be sufficiently high to be detected by GLAST. 



6 SUMMARY 

The current models of large scale AGN jets relate the non- 
thermal X-ray emission of distinct jet features, like the knots 
and hot spots, to the synchrotron or inverse Compton radi- 
ation of directly accelerated electrons. Both models however 
face serious problems. The inverse Compton or synchrotron- 
self Compton models typically fail on energetic grounds, 
unless one assumes that the X-ray emitting regions arc jet 
structures moving relativistically towards the observer. This 
is a very promising approach which seems to be applicable 
for the jets in 3C 273 and PKS 0637-752, but needs further 
inspections based on larger source statistics. 

While the weakness of the inverse Compton models orig- 
inates from the lack of sufficiently dense photon target fields, 
the electron synchrotron model has just an opposite prob- 
lem. It is "over-efficient" in the sense that the TeV electrons 
due to severe radiative losses have very short propagation 
lengths, and thus do not allow formation of diffuse X-ray 
emission on kpc scales, as it is observed by Chandra. A pos- 
sible solution could be that the electron acceleration takes 
place throughout entire volume of a knot or a hot spot. 
However the operation of huge, kpc size accelerators seems 
to be a non-trivial theoretical challenge. The secondary ori- 
gin of TeV electrons, produced homogeneously in the knots 
by relativistic protons due to interactions with the ambi- 
ent gas, seems an interesting possibility. But this hypothesis 
requires unacceptably large product of the gas density and 
the proton acceleration rate, unless we assume that the X- 
ray emitting regions are relativistically moving jet structures 
with Doppler factors S> 1. 

In the present work I propose an alternative mechanism 
for X-ray emission - the synchrotron radiation of extremely 
high energy protons with Bp > 10^^ eV accelerated in the 
in large-scale jet structures. The idea advanced in the pa- 
per is that it is possible to construct a realistic model which 
allows effective cooling of protons via synchrotron radiation 
on quite comfortable timescales of about 10^ — 10* yr, i.e. 
on timescales which provide effective propagation of pro- 
tons over the jet structures on kpc scales. This explains in a 
rather natural way not only the diffuse character of the ob- 
served X-ray emission, but also the broad range of spectral 
indices observed from different objects. Yet, the model pro- 
vides quite high radiation efficiencies if we allow relatively 
large magnetic fields in the knots and hot spots at the level 
1 mG. For example, the required proton acceleration rates 
ranges from lO"^ - lO"" crg/s in the 25" knot in 3C 120 and 
in the west hot spot of Pictor A, to 10"*^ — 10"*^ erg/s in the 
knots of powerful quasars 3C 273 and PKS 0637-752. For 
relativistic jots aligned with the line of sight these numbers 
can be reduced by two or three orders of magnitude. 

An essential fraction of the jet energy released in the 
form of extremely high energy protons eventually escapes 
the jot; the kpc scales of knots and hot spots are not suffi- 
cient for effective confinement of the most energetic particles 
with Bp > 10^^^ eV. The runaway protons interacting with 
the CMBR photons, as well as with the ambient gas, initiate 
non-negligible nonthermal radiation components in the clus- 
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ter environments of AGN and radiogalaxies. The signatures 
in the spatial and spectral distributions of this emission, 
which extends from radio wavelengths to very high energy 
7-rays, contain unique information about both the highest 
energy protons accelerated in the AGN jets and about the 
the intergalactic magnetic field in the extended intergalactic 
environments around powerful AGN and radiogalaxies. 
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